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The global first passage time density of a network is the probability that a random walker re- 
leased at a random site arrives at an absorbing trap at time T. We find simple expressions for the 
mean global first passage time (T) for five fractals: the d-dimensional Sierpinski gasket, T-fractal, 
hierarchical percolation model, Mandelbrot-Given curve and a deterministic tree. We also find an 
exact expression for the second moment (T 2 ) and show that the variance of the first passage time, 
Var(T), scales with the number of nodes within the fractal N such that Var(T) ~ N A / d , where d is 
the spectral dimension. 



1. INTRODUCTION 

The study of diffusion on statistically self-similar disor- 
dered media (e.g. percolation clusters and diffusion lim- 
ited aggregates) is simplified by modeling their structure 
by deterministic fractals [l[ . Although scaling theory has 
been central in studying diffusive processes on a variety 
of deterministic fractals @ , it has been noted that scal- 
ing laws do not provide a complete picture of dynamic 
phenomena on fractals [|[, and exact relationships are 
useful. 

Lately there has been interest in calculating the mean 
time (T) taken for a walker released from a randomly 
chosen node on a fractal to arrive at a trap. We call 
the underlying probability density of T the global first 
passage time (GFPT). It is found by a simple average 
of the node-to-trap first passage times over the entire 
fractal. The mean of the GFPT was first obtained for the 
Sierpinski gasket lattice by Kozak and Balakrishnan 0, 
Q , and an analogous study has been recently performed 
for the T-fractal p. The mean (T) has been associated 
Q with important questions such as the exploration of a 
random walker on a fractal lattice Q, providing further 
motivation for the study of global first passage times. 

In general, deterministic fractals have played a key role 
in explaining diffusion in disordered materials because 
their properties can be exactly studied. In this paper 
we use a renormalization method to derive an exact ex- 
pression for the full moment generating function of the 
GFPT. The method is exact, and can be applied to any 
finitely ramified deterministic fractal, making the explo- 
ration of higher moments of the GFPT possible. There 
is interest in higher moments of random walk quantities 
d, Q as such results give further insight into the nature 
of transport on fractals. Knowledge of these moments is 
also useful to investigate the variance of the bi-molecular 
reaction A + B — > A fi~rj||. 

The paper is organized as follows: In §2 we formulate 
the problem of continuum diffusion on networks of pipes 
as an algebraic system of equations involving Laplace 
transforms of the fluxes and concentrations at each node. 
We show how the equations are solved on a d-dimensional 
Sierpinski gasket lattice using a renormalization proce- 
dure. In §3, we derive the general formula for the Laplace 
transform of the GFPT density for a trap at an arbitrary 



node on the gasket. We demonstrate our method for 
the apex of the nth iteration of the d-dimcnsional gas- 
ket. Our results, for the mean of the GFPT, agree with 
prior work d, and are extended to show that, for any 
given trap, there is an exact analytic expression for the 
moment generating function of the GFPT. 

In §4, we demonstrate that the method can be used to 
calculate the GFPT moment generating function for at 
least one node on any structure provided that a renormal- 
ization procedure exists. A general formula is provided 
for the GFPT moment generating function for a finitely 
ramified deterministic fractal with two end nodes. As 
examples we consider the Mandelbrot-Given curve, de- 
terministic tree, hierarchical percolation model and the 
T-tree @ . The calculation of higher moments is demon- 
strated for the T-tree in 85. 



CONTINUUM DIFFUSION ON NETWORKS 



We begin by solving the diffusion equation 



(1) 



on a bar of length L with the non- homogeneous boundary 
conditions 

C(0,i) = Pi(t), C{L,t) = P 2 (t) 

and initial condition C (x, 0) = 0. Taking Laplace trans- 
forms of ((T|) and solving using the boundary and initial 
conditions we find 



, , sum (x-v/tT) sinh ((L — x) \/rT) 



sinh 



sinh (Ly/jj) 



(2) 

where wc denote Laplace transformed functions by lower 
case letters; 



c(x, s) = ££ (C(x, t)) 



C (x, t) e- st dt. 



The flux entering the bar at the point x = is F± = 
—D-S^C (x, t) and the flux entering the bar at x = L 

. For networks it is useful to 
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is F 2 - D-jj-C{x,t) 
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express these fluxes in terms of the concentrations pi at 
cither end. This gives rise to the equations 
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— sech ( 

J (3) 

where g = tanh (y^L) /V sD. We call the above a flux- 
concentration matrix. For simplicity we set D and L to 
be unity. 

Using these relations, diffusion on an arbitrary network 
can be formulated as a system of algebraic equations. To 
link two or more bars at a node we set the concentra- 
tions to be equal. As an example the system for the 1st 
iteration of the Sierpinski gasket lattice shown in Figure 
1 (a) is 
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Here ao = 1 an d bo = sech (^/s). 



(4) 



2.1. Initial and boundary conditions 

The fi in ((4]) are understood to represent the total flux 
entering node i. For mass conservation fi = 0, but if we 
take an instantaneous unit source at node i, then the flux 
entering the node is Fi = 5(t), where 5(t) is the Dirac 5 
function. This means that f t = Jzf (Fi) = jSf (£(*)) = 1. 
The matrix has a similar form to a simple finite differ- 
ence scheme for modeling diffusion/probability on a lat- 
tice. Prior studies [ll|, [l2|, G3 have also formulated the 
problem of diffusion on networks in terms of the Laplace 
transform of the concentration and flux at a node. Here, 
linear algebra is used to extend these ideas. 

It is useful to illustrate the formulation of two dif- 
ferent but closely related problems. If we take f T = 
(0,0,0,0,0,1) the solution p T = (p 1 ,p 2 ,p 3 ,P4,P5,P6) 
provides the Laplace transforms of the concentrations 
at the nodes if an instantaneous source is released at 
node 6, and mass is conserved at every node (including 
node 6). Next consider taking f T = (1, 0, 0, 0, 0, /g) and 
p T = (Pi,P2,P-},P4,P5, 0) and solving the 6 equations for 
p[ (i = 1, . . . , 5) and /g. In this case the equations corre- 
spond to an instantaneous source at node 1, a homoge- 
nous Dirichlet condition (i.e. a trap) at node 6, and mass 
conservation at nodes i = 1, . . . , 5. The quantity is 
just the flux entering at the "trap". It can be shown [14[ 
that this flux can be calculated from p\ and pe derived in 
the former problem, the result being /g = —pi/p^. Note 
the negative sign arises from our convention for defining 
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FIG. 1: (a) 1st iteration of the Sierpinski Gasket lattice, (b) 
0th iteration of the Sierpinski gasket lattice. The labeling 
used in the above (a) and (b) implies that the left corner of 
the nth iteration of the 2-d Sierpinski gasket lattice is denoted 
by 3 (3 n + 1) /2 - 2 n , the right corner by 3 (3™ + 1) /2 and the 
apex by 1. 



the flux. The result (a convolution in time) can also be 
derived using a continuum random walker argument. We 
use a generalization of this result in §3. 



2.2. First passage times and concentrations 

In order to analyze the flux-concentration matrix for 
higher order generations of the Sierpinski gasket we apply 
a renormalization method. We first show how the 6x6 
system (J4j) can be reduced to a 3 x 3 system associated 
with the 0th order Sierpinski gasket (See Figure 1 (b)) 





■fi- 




"2a -bo -bo' 




'Pi ' 


9 


h 




-bo 2a -bo 




P2 








. -bo -bo 2a . 




-P3 - 



(5) 



To do this we consider the three equations associated 
with nodes 2,3 and 5 in equation j4]) 
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where we have set f 2 — 0, fa = and /s = 0. Solving 
this under determined system gives 
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where po = f 2 -. Eliminating P2,P3 and p^ from ([4]) by 
using ([5]) gives 



FIRST PASSAGE TIME 
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Note that this flux-concentration matrix is of the same 
form as a flux-concentration matrix (|5|) of the triangle 
in Figure 1 (b). Continuing, we find that the system of 
equations for the nth iteration of the Sierpinski gasket 
has the form |(SJ) with ag and bo replaced by 



PrJ^n Pn 



b n -i (2p ra _i + 1) 
(4p„_i + l)(2 /9n _i- 



1 



(I) 



Here p n is the renormalization of the inverse FPT equa- 
tion on a 2-d Sierpinski gasket, which satisfies 

Pn = Pn-l (4/0„_i - 3) . 

A detailed discussion of the equation can be found in [l5|, 

MM- 

The recurrence relationship for the d-dimensional Sier- 
pinski gasket lattice is 



pnbn Pn 



where 



(dp„-i + 1) b n . 



{dpn-i - (d - 1)) (2dp n - X + (3 - d)) ' 



(8) 



p n (d) = 2d P 2 n _ 1 - 3 (d - 1) p„_a + (d - 2) (9) 

is the renormalization of the inverse FPT equation in d- 
dimensions. Equation © has been derived in [l9j 
however the above recurrence relationships and (jSJ) 
are new, and allow the concentration to be calculated 
everywhere within the nth iteration of the d-dimcnsional 
Sierpinski gasket lattice given that an instantaneous unit 
source is released at a corner of the structure. 

Before proceeding, a comment needs to be made about 
equation ©. The relation between the concentrations 
calculated in ([6]) holds for any iteration. By this it is 
understood that if P2,P3,P5 can be expressed in terms 
of pi,P4,pq through ([6]) using po for the 1st iteration of 
the Sierpinski gasket, then Pi,pe,pi3 can be expressed in 
terms of pi,pn, P15 through j6|) using p\ on the 2nd itera- 
tion of the 2-d Sierpinski gasket. In our notation (Figure 
1), node 13 is the mid base node equidistant to the nodes 
11 and 15, with the other nodes having a similar label- 
ing to Figure 1. The fact that the concentrations at the 
nodes can be expressed through (JS]) will be used when we 
derive a similar relation that holds for the concentrations 
at the nodes of other fractal structures in 54. 



Let Fkj (t) be the first passage time for a random walker 
released at node j to reach node k. The GFPT is given 
by 



F k {t) 



1 

N 



N+l 



where N + 1 is the total number of nodes and N the 
number of starting nodes on the fractal. Let Pkj, n be the 
concentration at a node k given that an instantaneous 
unit source is released at a node j, and there is no loss 
of mass at all nodes (the first problem discussed in §2.ip . 
The additional subscript n denotes the iteration of the 
fractal. As in t )2.1l wc consider a complementary problem 
where a trap is placed at node fc and a unit source is 
released at node j. The flux "entering" at k due to the 
source at j is 



fkj,' 



Pj k , n 
Pkk,n 



The Laplace transform of the GFPT for node k is 



N+l 



Zjj'=1 Pjk,n 



jk,n 



J'=l 



Pkk.n 



(10) 



The mean of the GFPT is 



tF k {t)dt 



dfk,n( S ) 



ds 



The form of fh, n 031 is simplified for the Sierpinski gas- 
ket through the use of conservation of probability. We 
integrate equation ((2|) for every pipe of the gasket and 
sum to find 



d {a k , n + 2 j3 k ,n) (cosh (yg - 1) 
Vssinh(^s) 



1 



(11) 



Here a k , n is the sum of concentrations at the corner 
nodes and the apex for the closed Sierpinski gasket 
and Pk,n is the sum of the concentrations at the other 
(d + l) ((d + 1)" - 1) /2 nodes. Note that (jTTJ is true be- 
cause there are d pipes stemming from the corner nodes 
and there are 2d pipes stemming from each interior node. 

Using the method of images we show that a k n = 
P\ k ,n-\- Here we have identified the apex of the nth iter- 
ation of the d-dimensional Sierpinski gasket by the label 
1. We consider the release of a walker at point k = 2 
of the 2nd iteration of the 2-d gasket (see Figure [2] (a)) 
to show the result, but the extension to arbitrary k and 
higher dimensions is clear. We first place an additional 5 
sources at the image points of k = 2 in the three lines of 
symmetry of the gasket (Figure [3 (b)). If we denote the 
sum of concentrations at the corner and apex nodes by 
a' 2 2 ) then it can easily be seen (through symmetry) that 



4 




FIG. 2: Illustration of the method of images argument used to prove the result a k ,n = Pik,n-i- The filled circles denote sources 
and the open circles are the concentrations of interest. 



&2,n = £*2 2/6- Now, by symmetry, there will be no flux 
between the three 1st iteration gaskets used to construct 
the 2nd iteration, and each problem (see Figure [3] (c)) 
can be treated alone. It is clear that the concentration 
at the apex of the reduced problem is 2pyi,\- Therefore, 
= 3(2p 12j i) and finally 0:2.2 = Pi2,i- For a source 
at a point which lies on a line of symmetry (such as a 
corner), only two copies of the source need to be consid- 
ered. In the above argument k is limited to values which 
appear in the n — 1th iteration. All other points in the 
nth iteration have an image point corresponding to an 
acceptable k value. 

Using (UT]) reduces the GFPT (10} to 
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where N (d) = ((d + 1)" +1 + d - lj /2 is the number of 
starting nodes. Setting ak, n — Pik.n-i m (fT2|) gives 

Pkk.n Sinh (Vs) 



fk 




Plk,n-1 



dy/s (cOSh (y/s) - 1) Pkk,: 



(13) 

As an example, consider a trap at the apex which re- 
quires evaluation of pu, n , which is provided in (|A3jl of 
Appendix [A] After simplification the moment generating 
function of the GFPT for a trap at the apex is 



fl.n — 



dp r , 



dp j + 1 



2N \ d Pn -(d-1) 



2 + U dPj - {d -i) 



(14) 



3.1. Discrete random walker 



Note that the diffusion equations and first passage 
times considered above correspond to a continuous- 



space/continuous-time random walker. We have chosen 
to formulate the problem in this way because it makes 
the working significantly simpler. For example, (|10p ap- 
plies to all networks, while it must be modified to ac- 
count for local coordination numbers in the discrete for- 
mulation. From fk,n{s) we can find results for discrete- 
time/discrete-space walkers using the simple transforma- 
tion discussed below. 

Let F^-(m) denote the probability that a walker re- 
leased at j arrives at k on the mth step for the first time, 
so the GFPT is given by 



N+l 



Ft 



N ^ ' 3 

j=hj^k 



F* k (m). 



Here, and below, wc denote quantities associated with a 
discrete walker by an asterisk. The moment generating 
function of F^ (m) is defined as 



77 1 A 



rn=0 



/y(A) is the counter part of fkj(s) descibed above (the 
iteration for the fractal is irrelevant here). The two quan- 
tities are related by 



fkjW = fkj ((arscch( e A )) : 



For example, on a single bar T^ 2 (m) = 5 m i , so /i 2 (A) = 
e A and for a continuum walker /i 2 (s) = l/po(s) = 
sech(y / s). A derivation of the transformation is given 
in Appendix [Bl 

Thus, the moment generating function of the GFPT 
for a discrete walker is given by 



/k(A) = fk ((arsech(e A )) : 



(15) 
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so we can write 



m=0 



r/A 



A=0 



Since the location of the trap is generally clear, we simply 
write the mean as (T*). To explicitly calculate F£ (m) , we 
form the probability generating function /^(ln(z)) and 
note that F£ (m) is the coefficient of the mth term of the 
Taylor series expansion of /£(ln(z)) about z = 0. 

For the remainder of the paper, we first calculate /fc(s), 
and then use (fT5|) to determine f£(\). For the Sierpinski 
gasket example, a Taylor series exansions of the moment 
generating function /^(A) ( (fl4|) in (fl5|) ) determines the 
moments of the GFPT: 



(d + 1)" ((d + 3) n+1 - l) + (d + 2) (d + 3)") d 2 



((d+l)" +1 + d-l) (d+2) 



A 

+0(A 2 ). 
(16) 



The coefficient of A is the mean GFPT which matches the 
results found in 0, [f| . The first and higher moments of 
the GFPT for alternative nodes can be calculated using 
our method. For example, the GFPT for a trap at the 
mid-base of the 2-d Sierpinski gasket is 



2N\ Pn ij2 ft -l Pn 



which allows us to calculate f% n (A) for the case of an end 
node. 

In the next section, we will always calculate fk,n(s) 
for node 1. The trick in determining pop is to express 
the concentrations at the interior nodes in terms of the 
concentrations at the exterior (end) nodes. This is rele- 
vant when calculating (JTUJ) in terms of the concentrations 
throughout the fractal, which can then be explicitly writ- 
ten as a function of the renormalization of the inverse 
FPT equation. /^„(A) can then be found from (T5]). As 
examples we calculate n (A) for four different fractals. 



(a) 



(b) 



FIG. 3: Deterministic tree (a) n = 1 (b) n = 0. The labeling 
of the nodes 1 to 5 remain the same regardless of the iteration. 



4.1. Deterministic Tree 

For the deterministic tree (Figure [3]) we calculate 
fk,n(s) for the end node 1. We can relate the concen- 
trations at the nodes 2 to 4, for the iteration n = 0, to 
the concentrations at nodes 1 and 5 in the following way: 



where k = 3 (3™" 1 + l) /2 - for n > 1, using the 

labeling in Figure [TJ For a discrete walker, we find 



2 (15™) - 3" + 11 (5"" 1 ) 
3 n+1 + 1 ' 



4. METHOD FOR OTHER RENORMALIZABLE 
STRUCTURES 

For completeness it is useful to determine whether 
fkn(^) f° r an y finitely ramified deterministic structure 
can be found. In the next section we are able to obtain 
at least one exact analytic expression. /^ n (A) for other 
structures can also be calculated depending on the coor- 
dination number of the nodes in the network. For the 
Sierpinski gasket a simple expression for f£ n (A) can be 
found for more than one node. This is due partly to the 
structure of the Sierpinski gasket, where every node is 
connected to 2d other nodes with the exception of the 
corners. In that case, conservation of probability simpli- 
fied the solution but such a formula might not exist for 
other structures. However, for all finitely ramified deter- 
ministic fractals there is always a renormalization scheme 
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Here p n +\ = 4/0 2 — 3 is the renormalization of the inverse 
FPT equation. We can deduce from (Ti~7|) that 
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and also that 



5 
i=l 



2p (l+p ) 
2p\-\ ' 



(18) 



(19) 



where Ki is the coordination number associated with the 
ith node. These equations are essential, as they can be 
used repeatedly on any iteration of the deterministic tree 
to express the concentrations at the interior nodes in 
terms of the concentrations at the nodes 1 and 5. This 
will now be demonstrated for iteration 1, shown in Fig- 
ure [31 which we label consistently with iteration (we 
are not concerned with the labeling of the new nodes). 
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We denote 02j to be the sum of the concentrations at the 
nodes between node 2 (i.e. the center node) and node j 
on iteration 1. Using (JT5|) 
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(Pll,l + P21,l) 



Repeating this process for all unlabeled nodes and sum- 
ming leads to 
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The final step is recalling (§2) that if the sum of concen- 
trations on the LHS of {T5J) gives the RHS of (JT5J) with p 
for iteration 0, then one can automatically assume that 
the LHS of(T5J) gives the RHS of (fill) with p being re- 
placed by pi for iteration 1. In doing so, one finds the 
sum of the concentrations at all non-labeled nodes on the 
fractal at iteration 1 is 



721 + #23 + #24 + #25 = 



(po + 2)(l + pi)pi 
(2pg - 1) (2pf - 1) 



Furthermore, one can show that all concentrations at any 
node at any iteration can be expressed in terms of the 
concentrations at nodes 1 and 5 with the sum of concen- 
trations at every node being 
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We know that there are exactly TV = 4 n+1 starting nodes, 
so /i, n (s) can be calculated from (| 10[) 



/i,»( s ) = ^ 
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Substituting for A ([To]) and taking Taylor series expan- 
sions gives 



3(2 ™-i )+ 6 8 i +n _17 



A + 0(A 2 ). 



More generally, one can construct a formula for any 
finitely ramified dctcrminstic fractal, provided that it has 
two end nodes each having a coordination number of 1. 
In this case, one can find the following quantities 



2/i (Po) = 



1 



Pn,o + Pi/i,o 



TABLE I: The first and second derivates of Pn(A), evalu- 
ated at A = 0, for the d- dimensional Sierpinski gasket, the 
Mandelbrot-Given curve (M-G), the T-tree, the deterministic 
tree (D) and the hierarchical percolation model (HP). 
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where a is an interior node, and 
V2 (po) 



PU,0 + Pi/1,0 ~[ 



N 



(21) 



For consistency we always let 1 and v represent the two 
end nodes of the nth iteration of the specific fractal. 
One can then show that the GFPT density for a ran- 
dom walker to arrive at node 1 can be expressed as 
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which becomes, using the fact that p v i,n/Pn,n = Pn+i> 

1 r . 

fl,n( S ) = 
n / n 
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(22) 



Here p n +i is the renormalization of the inverse FPT equa- 
tion of the nth iteration of a finitely ramified determinis- 
tic fractal with two end nodes 1 and v. To calculate (T*) 
for discrete walkers we use /o*(A)|a=o = T & n d the first 
derivative of p* (A) given in table I. The method will be 
demonstrated on three simple examples. 



4.2. T-tree 

Consider the T-tree shown in Figure[2J The decimation 
procedure assumes that flux leaves only through the end 
nodes 1 and 4 for every iteration n. The concentrations 
at nodes 2 and 3 for the iteration n = can be expressed 
in terms of the concentrations at nodes 1 and 4 in the 
following way: 
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FIG. 4: The T-tree (a) n = 1 (b) ra = 0. The labeling of the 
nodes remains the same for all iterations. 



Hence, the quantities defined in (fSU|) and (|2"Tj) are 

Po(l + Po) 



2/2 (Po) = 3po2/i (Po) = 3 



3/>o 2 -l 



We use the renormalization of the inverse FPT equation 
p n+1 = 3pl - 2 and N = 3" +1 in equation $22]), to find 



/l,n(s) 

n 3 



i 



Pn+1 



AT 

v 3p m 2 -l^ +1 " 3 Pi 2 -l 



(23) 



Substituting for A (fT5|) in (|23|) and taking Taylor series 
expansions gives 



-on+2 4 7 



A + 0(A 2 ). 



The mean coincides with the (T*) for a trap at the center 
(node 2) of the T-tree of iteration n— 1 (replace n by n— 1 
in the above formula), which agrees with [f|. 

For the next two examples we just state N, the renor- 
malization of the inverse FPT equation, and the two es- 
sential quantities yi(po),y2(po)- 



(a) (b) 

FIG. 5: Mandelbrot-Given curve (a) n = 1 (b) n = 0. 



3 



13/2 p^ 1 - 1/2 p n - 3 + 16 pr, 

2pl + 2p -l 



and 



2/i (A)) 



8p3-6p 2 -2p + l 

2/QQ (4 / o 2 +/7 -l) 
8^-6p 2 -2p + l' 



Substituting these expressions into equation 
in the usual way that 



, we find 



= ! + 4/7 1518(176")- 39(8 : ) + 187(22") A 



-0(A 2 ). 



The coefficient of A is the mean GFPT for a random 
walker to reach node 1 of the nth iteration of the 
Mandelbrot-Given curve. 



(b) 



(a) 



o 




FIG. 6: Hierarchicial percolation model (a) n = 1 (b) n = 0. 



4.4. Hierarchical percolation model 

The hierarchical percolation model shown in Figure [6] 
has the properties N = (2(4 n+1 ) + l)/3, p n+1 = p n {§p 2 n - 
7) /2 and 



V2(po) = (3po + l)2/i (Po) 



(3po + 1) 
3p -2 ' 



Using these results in equation (|22|) . transforming to A, 
and taking Taylor series expansions gives 



A* (A) = 1 + 2 



10 



n+l 



30(40 



n\ An 



1 + 2(4™+!) 



-A + 0(A 2 ). 



Here, as in all other examples, the coefficient of A is the 
mean of the GFPT for node 1 of the hierarchical perco- 
lation model. 



4.3. Mandelbrot-Given curve 



SCALING OF HIGHER ORDER MOMENTS 



The Mandelbrot-Given curve is shown in (Figure O . 
We find that N = (6 (8 n+1 ) + l) /7, p n+1 = -21 p n + 



To find higher moments of the GFPT, it is necessary to 
calculate higher derivatives of p* n (A). We have given the 
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relevant results in Table I. As an example, we calculate 
the second moment of the GFPT on a T-Tree. From the 
coefficient of A 2 in the Taylor expansion of /*(A) 
I15p we find that the second moment is given by 



((T* 



72 



12" 



16 
15 



2 n 



8928 
~L75~ 



36" 



208 
~25~ 



6" 



53 



175 

(24) 

It is interesting to examine how the first and second mo- 
ments, expressed as a function of the number of nodes 
Nt = N + 1, scale with the spectral dimension d [20I ]. 
Usually, d is determined through the analysis of the 
Green's function at the origin Go(t), which is known to 
scale as 

G (t) ~ t~*' 2 . 

In 0,11, 

the mean GFPT expressed in terms of the total 
number of nodes Nt on the Sierpinski gasket was re- 
lated to the spectral dimension d = 2 ln(d+ 1)/ \n(d + 3). 
This was done by re-arranging the formula for the num- 
ber of nodes on a rf-dimensional gasket, Nt = N + 1 = 
(d+1) ((d+l) n + 1) /2, to set (d+l) n = N T -^-l and 



(d+3) n = (N 7 



d+i 



T) 2 / d . Substituting these values into 



2/d 



(TlTJI and taking the limit for large N T gives (T*) ~ N, 

We perform a similar analysis for the second moment 
of the GFPT given in flU}. ' There are N T = 3" +1 + 1 
nodes available on the nth iteration of the T-fractal. By 

replacing 3" by {N T - 1) /3 and 2" by ((N T - 1) /3) 2/d_1 
in the expression for the second moment of the GFPT 
(EH) we find 



((T*) 2 ) = 6/5 (Nt - l) 4/3 - x (N T - 1) 2/J_1 + 

15 

248 , Ar ,,4/J 104 . . 2/d 53 

(N T - 1) 1 (N T - I) 1 H ■ 

175 v 1 ' 75 v 1 1 175 



(25) 



Asymptotically we find that ((T*f) ~ N T /d . Since 
we have an exact expression for the mean and second 
moment wc can show that the variance scales as 



Var(T*) 



N. 



4/d 



(26) 



All other structures considered above yield the same 
asymptotic behavior. 



6. CONCLUSION 

We have used the properties of the renormalization of 
the first passage time equation to find the first passage 
time moment generating function for a variety of fractals. 
From this, we were able to derive simple expressions for 
the first and second moments of the global first passage 
time. 

Our results can be applied to a diffusion controlled re- 
action experiment of the form A + B — > B. Here A is a 
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FIG. 7: The normalized probability for a random walker to 
arrive at a trap at node 1 on the mth jump for the 10th it- 
eration of the Sierpinski gasket (d — 2) (o), Sierpinski gasket 
(d = 3) (+) and the hierarchical percolation model (o). The 
values are calculated from the Taylor series of the GFPT gen- 
erating function f*(\n(z)). A symbolic algebra package was 
used. 



particle which performs a random walk on the structure 
until it is absorbed by a single immobile particle B. If 
there are an infinite number of A particles then the reac- 
tion rate at node k is the deterministic function Fj*(m). 
This reaction rate, at node k = 1, for the gasket in two 
and three dimensions, and the hierarchical percolation 
model, is shown for small time in Figure [7J The first and 
second moment we have calculated are the moments of 
the reaction rate. An alternate experiment would involve 
a single walker in which case the expected time of com- 
pletion is (T) with the outcome being cxtrcmly variable, 

the variance being given by Var(T*) ~ N^! d . 

Since the variance of the first passage time scales as 
([26f . we relate the variance to the mean by calculating 
the reduced moment R(N T ) = Va^T*) 1 / 2 / '(T*) . This 
type of approach has been used in the study of the vari- 
ance of the range of an n-step random walk 0, [2l| . As in 
[2~fl ] we expect that lirrijv T ->oo R(Nt) exists and remains 
finite, for every fractal structure with a spectral dimen- 
sion of d < 2. Using (fT4|) . wc find the reduced second 
moment for the Sierpinski gasket to be 

1/2 



lim R(N T ) = 1 + — 

Wr-00 v L > \ (d + 3)(d + A)d 

which decreases monotonically with d. 
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APPENDIX A: CORNER CONCENTRATIONS 

Here we show how the concentration at the apex and 
the corners of the nth iteration of the d-dimcnsional Sier- 
pinski gasket can be calculated if a walker is released at 
the apex. Wc define p c (i)i, n to be the concentration at the 
i = l..d corners of the gasket. Recall that the concentra- 
tion at the apex is given by pn >n . The flux-concentration 
matrix has the simple (circulant) form 







( da n -b n 
—b„ da n 



\ ( PllM \ 
Pc(l)l,n 



APPENDIX B: CONTINUUM TO DISCRETE 
TRANSFORMATIONS 

The results we derive for a continuous 
space/continuous time random walker can be ap- 
plied to a discrete walker by simple transformations. For 
a discrete walker the probability flux and probabilities 
at ends 1 and 2 of a bar are related by 

_ P?{m) Pj(m-l) 
l \ m ) — 

«1 K 2 



W \-b 



«1 K 2 

f?(nt) = gW- J ?("'- 1 

K 2 Kl 



■) 



b n da n j \ Pc(d)l,ra / 



This is the natural extension of equation ([5]) , with a n and 
b n defined in equation (JSj) . The solution of the equations 
are 



Pc(i)l, 



d (a n - b n ) {da„ + b n ) db n (p n - 1) (dp n + 1) 

(Al) 

for i = 1, . . . , d, and 



~fl ' 




1/ki 


—Z/K2 




Pi 


Jb. 




. -z/ki 


l/«2 . 




. Pz . 



Pll 



(da n - (d- l)b n ) 



(d Pn -(d-l)) 



d (o n - b n ) {da n + b n ) db n (Pn ~ 1) (dPn + 1) ' 

(A2) 

If we use the renormalization of first passage time equa- 
tion p n (d) = 2dp 2 n _ 1 - 3 (d - 1) Pn-i + (d - 2) (See ©) 
on the term (p n — 1) in pn. n (|A2|) . we find that 



where m > and re, is the coordination number of the 
ith node. These equations are quite formal and are easier 
to understand when they are combined for a network. 

Multiplying both equations by z n and summing from 
n = 0, . . . , co leads to 



(Bl) 



where Pi(z) = ^2^ =Q P* (jn) z m is the probability gen- 
erating function of P*(m). We have used the fact that 

define P*(— 1) = 0. Equation (|B1[) is the direct analog of 
©■ 

The transformations we need to derive are most easily 
shown by example. Using (|B1[) for the first iteration of 
the 2-d Sierpinski gasket gives; 



Pn, 



Using 



b n = 



we obtain 





(dp n -(d-l)) 




"A" 




" a - b f 


bo 
4 











db n (2dp n 


_i-d + 3) (p n _i - 1) (dp„ + 1) ' 




/a 






6o 
4 


bo 
2 


bo 
4 











/a 




bo b 
2 4 


«() 





bo 
4 


bo 
2 




(dp„_i + l)6„_i 




/4 




- b f 





«() 


bo 
4 





(dp n -l - 


(d-l))(2d /9n _ 1 + (3-d))' 




/5 




- b f 


bo 
4 


bn 
2 


"0 


bo 
2 








./e. 







fco 
4 





bo 
4 




(dp n - 


i-(d-l)) (dp n -(d-l)) 


where 


^0 = 


; 1 and 6o 




-z. 


We 


release 





Pi 




P2 




P3 




P4 




P5 




-P6 - 



Pll.n 

Note the similarities with equation (|A2p . If we repeat 
the above procedure we get the result 



Pll,n 



sinh (v 7 *) 



d^/s (cosh (y/s) — 1) 



nn 



d Pj -(d-l) 
dpj + 1 



(A3) 



The concentration at the corner points, defined in equa- 
tion (fAT|) . can be found using the same procedure 
used above. The result can be expressed as p c U)x,n = 
Pii,n/{dp n — (d — 1)). From the two preceeding results 
it can be shown that pii jn + dp c (;)i.n = Pii,n-i which 
agrees with the more general result given in 



(B2) 



at site 1 at t — and place a trap at site 6 with 
probability conservation at other nodes. This means 
P T = (Pi,P2,P3,P4,P5,0) and f T = (1,0,0,0,0,/e), 
where — /g is the probability generating function of the 
first passage time. The matrix equation given above can 
also be directly derived from the equations governing 
the probability of a discrete walker on the gasket (i.e. 
P x *(m + 1) = P 2 *M/ 4 + p 3*( m )/ 4 for n °de 1, etc). The 
initial and boundary conditions are represented in our 
flux vector. 

If we write the matrix equation (|B2[) as f = Ap and 
multiply by g — tanh(y / s)/y / s, we obtain 



gf = gAp. 



(B3) 
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Rcscaling p as 



2p[ 4p> 2 4p' 3 2p> 4 Ap' 5 



9 9 9 9 
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transforms the matrix equation (|B3|) to gi = A'p' . Be- 
cause of our rescaling of p. the matrix A' is now identical 
to the flux-concentration matrix defined in if z is re- 
placed by sech(y / s). Comparing the variables in both 
problems shows how first passage time and probability 
generating functions for discrete walkers are related to 
their continuum counterparts. 

More generally, for the quantities considered in this 
paper, the above example establishes the following trans- 
formations; 



1 tanh(^i) „ 
Pij{s) = 7= — Pij{z) 

Ki V s 



ch(v^) ' 



and fij(s) = fij(sech.(y/~s)). Recall / is the probability 
generating function of the first passage time. In order 
to relate / to the moment generating function /*, we 
make the further substitution /*(A) = /(e A ). Thus the 
quantities in the main text are given by 



fkjW = A,(e A ) = / fcj (arcsech(e A ) 2 ) 



which is the transformation used throughout the paper. 

We have not seen these transformations before, but 
they are analogous to the Montroll- Weiss transformations 
@, Page 254, Eq. 5.51]. which relate the properties of 
a continuous time discrete space random walker to those 
of a discrete time discrete space random walker. 
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